Statistical Isotropy violation of the CMB brightness fluctuations 
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Certain anomalies at large angular scales in the cosmic microwave background measured by 
WMAP have been suggested as possible evidence of breakdown of statistical isotropy (SI). SI vi- 
olation of cosmological perturbations is a generic feature of ultra large scale structure of the cosmos 
and breakdown of global symmetries. Most CMB photons free-stream to the present from the surface 
of last scattering. It is thus reasonable to expect statistical isotropy violation in the CMB photon 
distribution observed now to have originated from SI violation in the baryon-photon fluid at last 
scattering, in addition to anisotropy of the primordial power spectrum studied earlier in literature. 

We consider the generalized anisotropic brightness distribution fluctuations, A[k, n,r) (at con- 
formal time r) in contrast to the SI case where it is simply a function of \k\ and k ■ n. The bright- 
ness fluctuations expanded in Bipolar Spherical Harmonic (BipoSH) series, can then be written as 
Afff (k, t) where L > terms encode deviations from statistical isotropy. Violation of SI encoded 
in the present off-diagonal elements of the harmonic space correlation (ai m a\i m i) , equivalently, the 
BipoSH coefficients Afft 1 , are then related to the generalized BipoSH brightness fluctuation terms 
at present. We study the evolution of Af™ (k, r) from non-zero terms Af*£ (k, t s ) at last scattering, 
in the free streaming regime. We show that the terms with given BipoSH multipole, LM, evolve 
independently. Moreover, similar to the SI case, power at small spherical harmonic (SH) multipoles 
of Af^ 4 (k,T s ) at the last scattering, is transferred to Af^f 2 (k, r) at larger SH multipoles. The 
structural similarity is more apparent in the asymptotic expression for large values of the final SH 
multipoles. This formalism allows an elegant identification of any SI violation observed today to a 
possible origin in SI violating physics present in the baryon-photon fluid. This is illustrated for the 
known result of SI violating angular correlations due to the presence of a homogeneous magnetic 
field in the baryon-photon fluid. 



I. INTRODUCTION 

The Cosmic Microwave Background (CMB) anisotropy 
is a very powerful observational probe of cosmology. 
In standard cosmology, CMB anisotropy signal is ex- 
pected to be statistically isotropic, i.e., statistical expec- 
tation values of the temperature fluctuations AT(n) = 
Tli im a tmXim(n) are preserved under rotations of the sky. 
The condition for statistical isotropy (SI), in spherical 
harmonic space translates to a diagonal (ae m a}, m ,) — 
CiSu'Smm' where Ce is the widely used angular power 
spectrum of the CMB anisotropy. 

After the release of first year data of the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP), statistical isotropy 
of the CMB anisotropy attracted considerable attention. 
The study of full sky maps from the WMAP 5 year data 
and the very recent WMAP 7 year data 0], has 
led to some intriguing anomalies which seem to suggest 
that the assumption of statistical isotropy is broken on 
the largest angular scales Broken isotropy would 

have a profound implications for standard cosmological 
model as statistical isotropy underlies all cosmological 
inferences. 

It was pointed out that the suppression of power in the 
quadrupole and octopole are aligned in the form of the 



" axis of evil" [Tfj| - [l3 | . Further "multipole- vector" direc- 
tions associated with these multipoles (and some other 
low multipoles as well) appear to be anomalously corre- 
lated @ , [HI , [HI . There are indications of asymmetry in 
the power spectrum at low multipoles in opposite hemi- 
spheres, the "north-south asymmetry" Q, H, [T7rll9j . 
Possibly related, are the results of tests of Gaussianity 
that show asymmetry in the amplitude of the measured 
genus amplitude (at about 2 to 3 a significance) between 
the north and south galactic hemispheres [20|-|22j . Analy- 
sis of the distribution of extrema in WMAP sky maps has 
indicated non- Gaussianity, and to some extent, violation 

of sim. 

An observed map of CMB anisotropy, AT(n) obs con- 
tains the true CMB temperature AT(n) fluctuations, 
convolved with the beam and instrumental noise & 
foreground contaminations. Breakdown of statistical 
isotropy can occur in any of these parts and can be cat- 
egorized as 

• Theoretically motivated effects which are intrinsic 
to the true CMB sky, AT(n) include non-trivial 
cosmic topology 24 1, Bianchi models [25rf28j and 
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primordial magnetic fields [2!| , [3(| ■ A recent arti- 
cle [3l| claims that the solution to the cosmologi- 
cal vacuum energy can be explained as a result of 
the interaction of the infrared sector of the effec- 
tive theory of gravity with standard model fields. 
This theory predicts the violation of cosmological 
isotropy. 
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• Although a possible source of SI breakdown, resid- 
ual foreground contamination would need to be of 
order the intrinsic CMB temperature anisotropy to 
account for an appreciable effect [32M34] ]. 

• It would be erroneous to assume that the true CMB 
temperature fluctuations are completely extracted 
from the observed map. Observational artifacts 
such as non-circular beam, inhomogeneous noise 
correlation, residual striping patterns could be po- 
tential sources of SI breakdown. 

Violation of statistical isotropy of CMB anisotropy and 
its measurement has been discussed in literature earlier 
IHSl by defining an estimator where SI breakdown in 
an observed CMB anisotropy sky map is indicated by 
non zero value of this estimator. Studies have also been 
done by implementing a directional dependent inflation- 
ary power spectrum P{k) which gives rise to off-diagonal 
terms in the covariance matrix 41 1, [42] ]: spontaneous 



breakdown of SI in the CMB by a non-linear response to 
long-wavelength field fluctuations that appear as a gra- 
dient locally to the observer [43| or locally through a 
modulation field [44j]; incorporating an initial period of 
kinetic energy domination in single field inflation (45[. 

In this paper we present a new formulation that relates 
the breakdown of SI in the CMB photon fluctuations at 
last scattering, and evolving them to find the effect of 
the modes at present epoch hence the CMB (ae m a,£,i m >) 
today. We also find the Bipolar spherical harmonic coef- 
ficients (BipoSH) [3^, H^] which are linear combinations 
of off-diagonal elements of the covariance matrix. Bi- 
poSH expansion completely represents the information 
of the covariance matrix thus being the most general 
way of studying two point correlation functions of CMB 
anisotropy. These BipoSH coefficients are mathemati- 
cally complete measures of SI violation on a sphere. 



II. REVIEW OF STATISTICALLY ISOTROPIC 
CMB BRIGHTNESS FLUCTUATIONS 

A. Boltzmann equations, inhomogeneities and 
anisotropies 

In the smooth background universe, thermalized pho- 
tons being distributed homogeneously and isotropically, 
the temperature T is independent of x and direction of 
propagation p respectively. To describe perturbations 
about this smooth universe, we allow inhomogeneities in 
the photon distribution and anisotropies. 

Before recombination, z roc w 1100, the photons were 
tightly coupled to the electrons and protons; all together 
they can be described as a single fluid, the baryon-photon 
fluid. After recombination, photons free-stream from the 
surface of last scattering to the present epoch. 

Given the cosmological perturbations to the photons 
at recombination, one can predict the anisotropy spec- 
trum today. The main motivation next is to relate the 



moments today to the moments at recombination using 
the photon distribution function. 



B. Fluctuations of CMB photon distribution 

In the Boltzmann equation for photons df/dt = C[f], 
we expand the photon distribution function f(x,p, n, r) 
about its zero-order Bose-Einstein value T(r) [46| where 
p=pn. The distribution function of the photons changes 
with the perturbed temperature as 



f(x,p,h,r) 



exp 



T(r)[l + A(x,ri,r)] 



The perturbation to the distribution function is charac- 
terized by A = ST/T termed as CMB brightness fluc- 
tuations henceforth. Since the perturbation A is small, 
we can expand f(x,p, h, r) keeping only terms up to first 
order to get 



A(a ? , h, t) 



df° 
dlrip 



Sf, 



(2) 



where / is the zero-order photon distribution function. 
A(x, ft, t) depends on x, h and r and not on the mag- 
nitude of momentum p\ this is a valid assumption since 
the temperature of the plasma is very small compared to 
the rest energy of the electrons which undergo scatter- 
ing, elastic Thomson scattering has negligible effect on 
the magnitude of the photon momentum. 

Perturbations to the CMB remain small at all cosmolog- 
ical epochs; evolution of the largest scales being in the 
linear regime. In solving the linear evolution equations, 
it is simplest to work with Fourier transforms since every 
Fourier mode evolves independently. 



A(x, n, t) 



d 3 k 
(2tt) 3 



Ak-x 



A(k,h,r)4>(k) 



(3) 



where cj>(k) is the primordial density fluctuations. 
With statistical isotropy assumption 

A(fc,n,r) = A(k,k-n,r) 

= ]TH) £ (2£+l)A £ (fc,T)P,(fc.n) 

1 

= 47r^(-i)^A^(fc,r)^ m (fc)^ m (n). (4) 
im 

Ag(k,r) are the moments of the CMB brightness fluctu- 
ation. The monopole Ao is related to the density per- 
turbations while the dipole Ai oc h ■ v, gives the velocity 
term for baryons. 
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C. Correlations 

The observed anisotropy in multipole space can be 
written in terms of the CMB brightness fluctuation as 

A(x = 0, ft, t) = ^2a em Y em (n) . (5) 

Using the orthonormality property of the spherical har- 
monics Yi m , the SH coefficients ai m become 

/ii 3 k -» ~ 
_^(fc)A,(fc,r)y/ m (fc). (6) 

The angular correlation can be expressed as 

/ * \ a f & Pb (^) I A / , \ I 9 f c 

{at m a e , m ,} = 47r / — — — \At(k,T)\ Owb mm > 

/dk 
T P (fc)|A,(fc,r)| 2 , (7) 

where correlation of the primordial density fluctuations 
{cj)(k)cj)*(k')) = P (k)S(k - k') and V (k) = k 3 P (k)/2n 2 
is the primordial power spectrum per logarithmic interval 
generated by inflationary model and the second term is 
the radiative transport kernel in the post-recombination 
universe given by cosmological parameters. 

D. Evolution of CMB brightness fluctuation in the 
free-streaming regime 

The evolution of A(k,ft, r) in the free-streaming 
regime can be written as 

A(fc, ft, t) = e l£ -"( T - T =) A(fc, ft, t s ), (8) 

where r is well inside the free-streaming regime i.e. 
r s < r < To, t s and To being the conformal time at last 
scattering and today respectively Ho- 
using the expansion 

e ik.fiA T = J2(-i) l (2l + l)i,(fcAr)fl(* • ft), (9) 
i 

and defining At = r — t s , the evolution equation for 
Ai(k,r) can be written as 

Ae(k,r) = ^2(-iy+ i '- l (2£' + l)j l (kAr) 

w 

x [C l ° e , ] 2 A e ,(k, Ts ). (10) 

ji(kAr) and Pi(k ■ ft) are the £ th order spherical Bessel 
function and Legendre Polynomial, respectively. Equa- 
tion (fTUf is the well known "free streaming " equa- 
tion in CMB literature [47] . C^^ 2JTi2 are the Clebsch- 
Gordan coefficient which satisfies the triangle inequalities 
([HI) putting a constraint \£\ — £2] < L < £\ + £2 and 
mi +m2 — M. 



III. STATISTICAL ISOTROPY BREAKDOWN 
IN THE CMB BRIGHTNESS FLUCTUATION 

In this paper we take into account the SI violation of 
the CMB anisotropy which is seeded due to the inher- 
ent SI breakdown in the CMB photon distribution. We 
consider the general form of the CMB brightness fluc- 
tuation, allowing for anisotropy in k, i.e., A(fc, ft,r) ^ 
A(fc, k ■ ft, t). 



A. Generalized CMB brightness fluctuations 

The most general CMB brightness fluctuation is not 
simply a function of |fc| and k ■ ft. In this case the 
physical situation of anisotropic fluctuations demand the 
brightness fluctuations to be expanded in Bipolar Spher- 
ical Harmonic series (not just a Legendre series as in the 
statistical isotropic case). The brightness fluctuation in 
multipole space is Af_^/ 2 (k, t) where L > term incorpo- 
rates deviation from statistical isotropy. 

A(k,n,r) = 4ir ^ (i)^ yj £ x + £ 2 + 1 

£ 1 i 2 LM 

x AfX (*, t) {Y tl (k) ® Y l2 (ft)} lm 
= 4ir J2 Pt^^%(k,T) 

£ 1 l 2 LM 

x {Y tl (k)®Y^(ft)} LM , (11) 

where fii x £ 2 = (i) 5 \Jt\ + 12 + 1 has been defined for 

convenient notational simplicity. 

The tensor product in BipoSH function is defined as 

{Y ei (k)®Yt 2 (n)} LM = C trLi. m2 

77117712 

x Y eimi (k)Yi 

2 m 2 (ft)- 

The pre-factors in equation (llip are the normalization 
terms associated with the CMB brightness fluctuation. 
The deviations from statistical isotropy is associated with 
non-zero values of L, L > 0. To check for the statistical 
isotropy limit i.e. L = & M = 0, we use equation 
(8.5.1) from 

to recover equation dU in section (IIIBI) . 



B. Angular Correlations 

Starting with equation ([5]) and the Fourier Transform 
relation from equation ([3]), the SH coefficients ai m for the 
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general case are 

* ' iimiLM 

The angular correlations turn out to be 



{ae m a g , m , 



(4tt 



(2tt) 3 (2tt) 3 

E E faitPht' 



dmiLM l 2 m 2 L'M' 
X ^e 1 m 1 em^e 2 m 2 i'm' ^iroi (^)^ 2T n 2 

x A^(fc,r)[Ai;'f' (fc'.r)]*. (14) 

The most general power spectrum (under statistical ho- 
mogeneity) depends on the direction fc, 

(^(fe)0*(fc')) = P(k)6(k-k>). (15) 

Further, it is useful to parametrize the directional depen- 
dence of fc in P(k) as |42J 

i 

l>0m=-l 



P(fc) = P (k) 



(16) 



where the first term with [ = represents the statistical 
homogeneous and isotropic primordial power spectrum. 
For a directional dependent power spectrum, the angular 
correlations of temperature anisotropy can be written as 



{ai m a e r. 



= 47T 



k 2 dk 



2tt 2 

x [A0*'(fc,T>] 
\LMi 



Po(k) J2 

IimiLM L' M' 

E ° 2 9tm(k) 



DiA^(k,r 



x ^{k,T)[A£f{k,T)]*yli™tM >( 17 ) 



where 



a a* pLM pL'M' 
Pill Ptit' ^l imi lm WimiC' 

-iLM 



y\ Q £>* pl^lVl pL M 

t->i — Pi x t Pf, 2 £> ^e 1 m 1 e m ^i 2 m 2 i'- 



yl 2 m 2 



dQ- k Y iimi (k)Y; 2m2 (k)Y lm (k) 

n fl[ 



47T Hi. 



pt 2 pf, 2 m 2 
"^lOttr^iroilm' 



where we have used the expression for integral of three 
spherical harmonics as in equation (5.9.4) from (48[. 
Here U h e 2 ...e n = [(2£i + l)(2£ 2 + l)-(2i n + l)]i has been 
defined for convenient notational simplicity. 
The case L — reduces equation (fT7|) to that in the 
analysis [42j 



App,SwS r . 



+ (-l/' +m ']>>«' C 



(18) 



where statistical anisotropy is quantified by the BipoSH 
coefficients [36|, [38l - l40| , defined as a tensor product of the 
spherical harmonic coefficients ag m and ai> m r, 

MP — E ( aema i'm,') ("I)™ Qml'-m' 



JN 



= {a e <g>ai<} 
Aft = (-l) e U e C e 



(19) 



Here, Ci is the usual CMB power spectrum for the SI 
case. Directional dependent P(k) [42|], introduces the 
second term in equation (fT8f 



Aim 



Tl e U e 

n, 



[m 



00 



(20) 



where a corresponding tensor product in Bipolar har- 
monic space for the indices L and L' of CMB brightness 
fluctuations is defined as 



JN 

LML'-M' 



MM' 

A^(fc,r)[A^'(/ C> r)]*.(21) 



In general for SI violations (L > 0), the BipoSH coeffi- 
cients can be expressed using the angular correlations in 
equation (fT7|) as shown in Appendix [S] as 



.4 



JN 



k 2 dk 



Po(k) 



E n^^-i^'+^-^w? 



f 2 c 



J h I 



£ 3 m 3 I ^ ^2 



(22) 



l.m 



The first term on the right-hand side of equation (|22[) . 
is the contribution due to statistically isotropic primor- 
dial power spectrum i.e. [ = terms in equation (|16l) . 
The second term gives the contribution to the Bipo- 
lar coefficients due to [ > in equation (|16l) . The 
term in the first braces is the Wigner-9j symbol [48[ 
which is related to the coefficients of transformations 
between different coupling schemes of four angular mo- 
menta and satisfies the triangular conditions for the tri- 
ads (J£ 3 l), (tUx), (£'L'£ 2 ), (Jit), (£ 3 LL') and (U^)- 
The SI case L — reduces the BipoSH coefficients in 
equation (p?2")) to equation (|2"0) . 

To evaluate the first term (A'^^j ( in equation (f2"2"j). we 
consider statistically isotropic primordial perturbations 
(4>(k)(f>* (k')) — Po(k)5(k — k') and express the angular 
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correlations as 



IV. EVOLUTION IN THE FREE-STREAMING 
REGIME 



(ae m a}, m ,) = in 



k 2 dk 



2tt 2 



Po(k) 



E foil' 

ImnLML'M' 
LM/i _\ rXL'M', 



x A^(fc,r)[A^ (k,r)Y 

nLM nL'M' 
x ^imitoiWimifm'- 



A. Generalized evolution equation 

We find the moments of the CMB brightness fluctua- 
tion to be 



(23) 



1 1 



dflfi dQt A(k, n, t) 



To check for the statistical isotropic case, we put L = 
& M = and recover equation in section piCp . 

As shown in detail in Appendix [A] the BipoSH coeffi- 
cients in equation (j!9[) . can be expressed using the angu- 
lar correlations in equation (|23p for statistically isotropic 
primordial perturbations as 



4tt/^ 2 J J ~" " ^ k 

x {r fl (fc)®r, 2 (n)} LM , (27) 



starting with equation (jlll) and the orthonormality con- 
dition of BipoSH, 



dn n dn k {Y tl {k)®Y t M}LM{YiM®Y ti {h)Y LIM , 



5e,e, Se e, Sr L ,S 



4 OLL'OMM* 



(28) 



(Me?) i=o 



t+L+L' + J 



x {A^SA^,} (24) 



The first term in braces is the Wigner-6j symbol which is 
related to the coefficients of transformations between dif- 
ferent coupling schemes of three angular momenta. These 
vanish unless the triangular conditions |48| are fulfilled 
for the triads (LJL 1 ), (L't'l x ), (J£'£) and {£Li x ). 

We consider low Bipolar deviations from SI i.e. L, L' and 
J « £,£\£i in equation (|24|) and use the asymptotic 
relation for Wigner-6j functions given by equation (9.9.1) 



In the free streaming regime, using the plane wave ap- 
proximation as in equations (|8|) and ([9]), the most general 
evolution equation turns out to be 

r) = ^— E H)' n« j t {kAr) p i3ii 

P tlt2 Ua£iL'M' 

KLMr, T \ C LM nL'M' 

mim2m3m4 

x jj dQ n dk P e (k ■ h) y/ imi (k) y; ama (n) 
x Y l3ma (k)Y u 

)■ (29) 

As shown in detail in Appendix [B] this can be further 
simplified to 



from (48|. We find that the asymptotic limit to these A^,(fc,r) = ^ (-i) 1 Heu^ 



Ml I -,\l 3 +U+L 



BipoSH coefficients are 



(-1) 



(Me )i=o 



47r / ^ Po(fc) E ( - x 

J i\LL' 



ii+L+L' 



x h(kAr)AfX (k, t s ) C&° l0 Q&o 
. .£iL£ 2 



- Ull ' r r* r J( - e ~ e "> 



x {A^®Ai;V} jjv 



(30) 



(25) 



The generalized evolution equation thus can be expressed 
so as to structurally resemble the evolution equation for 
the SI case in equation (|T0l) . 



For diagonal brightness fluctuations i.e £\ = £ and £\ = 
£', the BipoSH coefficients themselves turn out to be di- 
agonal in multipole space 



Kt( k ^) = E C{L 1 £,£ l ,£ 2 ) jt(kAr) 



x A%t(k,r t 



(31) 



where 



(4« )l=o 



J T.T.I J 



LL 

x C'loL'O \j^U ® ^« } 



JiV 



(26) 



feO ^>£ 2 J ^1 L £ 2 1 



v ^«3U /i 



(32) 
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(a) Case 1: Plot of C(L,£,£ 1 ,£ 2 ) for &fi M (Jfc, r s ); with initial 
SH multipoles £3 = £4 = 1, the final SH multipoles being 
£ = £1 — 1(= £2 — 1) and the index of deviation from statistical 
isotropy L = (blue, solid), L = 1 (red, dashed) and L = 2 
(black, dot-dashed) respectively. 
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(b) Case 2: Plot of C(L, £, £1 , £2) for A% 2 M (k,r s ); with initial 
multipoles £3 = £4 = 2, the final multipoles being £ = £i(= £2) 
and the index of deviation from statistical isotropy L = 
(blue, solid), L = 1 (red, dashed) and L = 2 (black, 
dot-dashed) respectively. 



FIG. 1: Evolution of the coefficient of the spherical Bessel 
functions, C(L, I, li, £2) with multipole moment £1 as in equa- 
tion (|31[) . for non-zero diagonal terms of a unit normalized 

Af 3 f 4 (fc,r a ) 



Setting L — 0, M — we recover the statistical isotropic 
case as in equation (jTUl) . section (|IID|) . 
The transfer of power of the statistical anisotropic terms 
to higher SH multipoles £\ & £2 due to free streaming 
is illustrated in figures ([I]) and ([2]). Starting from a unit 
normalized A|^ 4 (k,r s ), we plot the evolution of the coef- 
ficients in equation (p?2"j). C(L, £, £\ 1 £2) with ^1 for specific 
values of the SH multipole moments £3 and £4. We find 
that the values of £ and ^2 ar e constrained by the values 
of £\ due to the triangular inequalities of the Clebsch- 
Gordon coefficients. 

Figure Q] shows the evolution of these coefficients for 
diagonal, unit normalized CMB brightness fluctuations 
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(a) Case 1: Plot of C{L,£,£ 1 ,£ 2 ) for Af M (fc,r s ); with initial 

multipoles £3 = 1, £4 = 0, the final multipoles being 
£ = £ 2 = £,_ + 1 (blue, solid) and £ = £ 2 = £1 - 1 (red, dashed) 

respectively. The terms are non-zero only when the index of 
deviation from statistical isotropy L = \£g — £i \ = 1. 




2 4 6 8 10 

h 

(b) Case 2: Plot of C{L,l,li,l 2 ) for A^ M (fc,T s ); with initial 
multipoles £3 = 2, £4 = 0, the final multipoles being 
e = £ 2 =e 1 (blue, solid) and £ = £2 = £1 - 2 (red, dashed) 
respectively. The terms are non-zero only when the index of 
deviation from statistical isotropy L = \£^ — £^\ = 2. 



FIG. 2: Evolution of the coefficient of the spherical Bessel 
functions, C(L, I, £1, £2) with multipole moment £\ as in equa- 
tion (I31[) . for non-zero off-diagonal terms of a unit normalized 

T.) 

namely Af' 1 M (A:,T s ) and Af 2 M (fc,r s ) with L = 0,1,2, for 
possible values of £ and £2 ■ 

Figure [2] shows the evolution of the coefficients for off- 
diagonal, unit normalized Af' M (A:,T s ) and (k,T s ), 
for possible values of L,£ and £2- The coefficients for off- 
diagonal A^ 3 A | 4 (/c, t s ) vanish for L = i.e. the statistical 
isotropic case. 



B. SI violation at large multipoles 

The ability to measure violation of statistical isotropy 
at low SH multipoles is largely compromised by cosmic 
variance. At larger SH multipoles, the effect due to viola- 
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tion of SI would become prominent. The previous section 
shows that power is transferred from small to large SH 
multipoles during free-streaming. This opens the door 
to more readily measurable SI violations arising from SI 
violation induced due to physical processes (e.g., pres- 
ence of magnetic fields, or, other breakdown of rotational 
symmetries) in the baryon-photon plasma. 

Due to the tight coupling in the baryon-photon plasma 
prior to r s , primordial SI violation can be expected to be 
limited to small SH multipoles. It is illuminating then 
to obtain an expression for the free-streaming of BipoSH 
brightness fluctuations at small SH multipole moments at 
the last scattering to large SH multipoles at the present 
epoch for which we essentially evaluate the asymptotic 
limit of the CMB brightness fluctuations today. 
In equation (l3~Tj) we take £\ & 1% to be large compared to 
£3 &z £4. From symmetry of Wigner-6j symbols 



(33) 



We evaluate the asymptotic limit with arbitrary values 
of £, £1, £2, £3 and £4 using equation (9.9.1) from [48[ as, 

£4 L £3 \ (-l)^+L+d+c 

R+dR+eR+f 





t 


M 


/ 4 


L 


£,} 


1 £4 


L 


M 




£ 


£2} 



x C 



2£U t3 

ts{d-e) 

UU-e)L{d-f) 



(34) 



where R is large and 



d = ti-R 
e = e-R 

f = h-R- 

Thus for £\,£i » £3, £4, equation ([34]) becomes 

£4 L £3 1 . - « 



■4 ^ «3 

£1 £ £ 2 



21 IL 



-C 



Putting the expression for the Wigner-6j symbols from 
equation (|35|) . the asymptotic limit of the generalized 
evolution equation ([3Tj) takes the form 



fri K V2I: 



(-if 



x j e (kAT)Af^f 4 (k, t s ) C^° i0 C^j° 40 

x C Ull ~ l) (36) 

Equation (|36[) depicts how power in SI violating terms 
at small SH multipoles £3, £4 at the last scattering, free 
stream to higher SH multipoles £\ , £2 at the present 
epoch. Note the structural similarity to equation (fTUj) 
for statistical isotropic case. 

It is useful to provide explicit expressions for equation 
(1361) in two particular cases, when the asymptotic mo- 
ments of the CMB brightness fluctuation contains only 
diagonal terms and off-diagonal terms respectively. 



The evolution equation which involves only the diagonal 
terms of the moments of the brightness fluctuation at last 
scattering are 



A£l(fc,r) = ^2 C\ (L, £, tx) je(kAr) 

£ 



(37) 



with 

CiOM,4) 



Tip 



21 



(2€)!(2& + l) 



{2£i 



(24)! 



(24) ! 



(2£ 3 -L)l{2£ 3 +L+l)\ 



(38) 



where £,£1 » \£\ — £\. The details are given in Ap- 
pendix [C] The term under the square-root captures the 
L dependence of the free-streaming of L > terms. 

The evolution equation which involves only the off- 
diagonal terms of the moments of the brightness fluc- 
tuation at last scattering are 



x A^_ t ^_ t %T s )8 LM ^ (39) 



with 

C 2 (£,£i,£ 2 ) 



2i 



h+£ 2 -2£+l 
h+£ 2 + l 



(2£)\£ l \£ 2 \ 



(£!) 2 (£ 1 -£)\(£2-£.y. 
[2£ x - 2£+l)\{2£ 2 - 2£+l)\ 



(24 + l)!(2£ 2 + l)! 



.(40) 



where £, £ n » \£ n - £\, L with n = 1,2. The details 
are given in Appendix |D| As in equation (|38p . the term 
under the square-root captures the L dependence of the 
free-streaming of L > terms. 



C. SI violating physical effects at last scattering 

The patterns of the CMB temperature field i.e. the an- 
gular correlations observed today are traced back to in- 
homogeneities at the last scattering surface. In the tight 
coupling regime of the baryon-photon fluid, one expects 
power only at small SH multipoles of the Af^f (k, r s ). 
The generalized evolution equation of the CMB bright- 
ness fluctuations (f2"9"| free-streams this power at small SH 
multipoles in both SI and non-SI moments to correspond- 
ing SI and non-SI moments with same Bipolar moment 
L. Any observed violation of SI today is easier to inter- 
pret as generalized moments arising due to simple physics 
just beyond the fluid approximation regime. In this sec- 
tion, we illustrate this point explicitly for SI violation in 
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the CMB anisotropy in the presence of a homogeneous 
magnetic field at last scattering. 

The SH coefficients a^ m can be expressed in terms of 
the CMB brightness fluctuations at last scattering using 
equation (IT31 and the generalized evolution equation (j3Tj) 
as 

aim = 47T / <t>(k) 22 ^ C txm x im Y tirm{k) 

x Y, C(L,£ 2 ,£ X ,£) j, 2 (kAr) AfX(k,r s ), (41) 

where C(L,£2,£i,£) is defined in equation (|32|) . 
Hence, in general (a,£m a t'm') correlations measured at 
present are related to (AfY(k, T a )A| ' Jf (k,T s )} correla- 
tion between the generalized Boltzmann fluctuations at 
the last scattering surface. 

In particular, SI violation encoded in the off-diagonal cor- 
relation {at m ai' m t) (and non-zero BipoSH Af^ 1 , L > 0) 
are related as 

A$ ~ {A^(fc, r.) ® Ag e {k, t s )} jn (42) 

as in equations (f24l [25l [26]) or more generally for different 
Bipolar coefficients J' and N' as in equation (|2"2")l when 
the power spectrum is also anisotropic. 
Using equation (fTTj) , correlations of the CMB brightness 
fluctuations are 

(A(k, n, r s )A(k, n>, r s )) = (4n) 2 ]T £ ft 3 , 4 

x A^(fc,r s )A0/'(fc,r s ) 

x {Y tl (jt)®Y ta {h)} 

x \Y e3 (k)(g>Y e Jn')\ . (43) 

We illustrate the generality and power of our formalism 
using the case for a uniform magnetic field. We show the 
correlations of the CMB brightness fluctuations in this 
particular case is sourced by the bipolar dipole [L = 1) 
terms of equation (fTTj) with £\ = £2 = 1 where 

A{k,n,T s ) = 47riV3y)A^(k,T a ){H(fc)®ri(ft)} 

' — ' L J 1M 

M 

= 3i \/|E A ^( fc ' T «) (44) 

M 

with M = { — 1,0, +1}. Here (k X n)M is the usual cross- 
product written as irreducible products of the rotation 
group [HJ|. 

Using standard vector identity [48| and equation © 

(n • n')(fc • fc) — (n • k)(n' ■ k) 

= (hx k) ■ (n' X k), (45) 



the temperature correlations in presence of a uniform 
magnetic field as discussed in [29[ (see equation A3) are 
given by 

(A(fc,n,r s )A(fc,n',r s )) oc A}f(fc,r s ) A{f (fc,r s ) 

x (h X k) ■ (ri X k). (46) 

Here the bipolar dipole terms of the CMB brightness fluc- 
tuation A\^(k,T s ) Alf 1 (k,r s ) in equation (j44|) . encap- 
sulates the source term due to the presence of a uniform 
magnetic field [29(. 

In equation (|41"j) . the bipolar dipole terms A\f { (k,r s ) 
gives rise to SH coefficients 

r (pk - * 

Ofm = 47T / — — 3 (j)(k) 22 Pll C lmUm Y lmA k ) 

x ^[(...)^-i(fcAr) + (...)im(AAr)] 

x A}f(fc,r s ). (47) 
The angular correlations in this case turn out to be 

(ae m a}, m ,) = (•■•) ^<^mm' + (■••) ^'<±2<Smm'' (48) 

It is interesting to note that the known diagonal (£! = 
£) and off-diagonal {£' = £ ± 2) correlations in presence 
of a homogeneous magnetic field [29|, |30j , can be easily 
recovered in our approach. Work is in progress to relate 
other cases of SI violation originating in the physics at 
the last scattering surface using this formalism. 

V. CONCLUSIONS 

The search for subtle statistical isotropy breakdown in 
the universe is highly motivated by numerous theoreti- 
cal scenarios. The fluctuations in the cosmic microwave 
background is the arguably the most promising observa- 
tional probe of the SI of the universe. The violation of 
SI could have its origin not only in in anisotropic primor- 
dial power spectrum, but also in the SI violation in the 
fluctuations of the baryon-photon fluid at last scattering. 
SI deviations generated by a general form of anisotropic 
primordial power spectrum for isotropic Boltzmann func- 
tions has been studied in the recent literature [42|]. This 
paper includes this equally important possibility of a gen- 
eral SI breakdown in the CMB photon distribution func- 
tion. We study the generalized case of SI violation in 
terms of Bipolar spherical harmonic (BipoSH) bright- 
ness fluctuations, substantially extending the scope of 
origin of SI violation solely from the anisotropic primor- 
dial power spectrum. 

The breakdown of SI in the CMB brightness fluctu- 
ation results in off-diagonal terms in the SH space an- 
gular correlations (ae m ae> m ')-> or > equivalently, in the co- 
efficients the Bipolar spherical harmonic (BipoSH) rep- 
resentation [3|1 |39| . We relate the measurable BipoSH 
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coefficients to SI deviations in the baryon-photon fluid, as 
well as, the primordial power spectrum. The observable 
BipoSH coefficients can be compactly expressed in terms 
BipoSH brightness fluctuations terms through products 
of standard Clebsch-Gordon coefficients and a Wigner-9j 
function. We also present the expression for the simpler 
case of an isotropic primordial power spectrum, where the 
BipoSH coefficients turn out to be given through a com- 
pact combination of a Wigner-6j symbol and a Clebsch- 
Gordon coefficient. We also provide the large SH mul- 
tipole limit for these coefficients for the terms encoding 
deviations from SI at low BipoSH multipoles. 

We obtain the generalized free-stream evolution equa- 
tion for the SI violation encoded in terms of the Bi- 
poSH brightness fluctuations introduced in our work. 
We demonstrate that different modes BipoSH brightness 
fluctuations at the present epoch have to evolve from 
same Bipolar modes at the last scattering. The moments 
of the CMB brightness fluctuations Af¥(k, r s ) at last 
scattering are expected to non-zero at small values of SH 
multipoles l\ and £2 due to tight coupling. However, 
our results show that the power in these SI violating 
terms at low SH multipoles would be transferred dur- 
ing free-stream evolution to higher multipoles £3 and £^ 
in (fc, t) at the present epoch. This is akin to the 
well known free-streaming evolution of power in the SI 
brightness fluctuation at low SH multipole power at last 
scattering to large SH multipole at present. For clearly 
highlighting the structural similarity, we present the evo- 
lution of BipoSH brightness fluctuations in the asymp- 
totic case of large values of the final SH multipoles today 
relative to the initial SH multipoles at last scattering. 

While many of the claimed observational evidence of 
SI breakdown, such as the "axis of evil", "north-south 
asymmetry" etc., pertain to relatively small values of the 



SH multipoles where the significance is largely obscured 
by dominance of cosmic variance. However, SI violation 
at small SH multipoles in the baryon-photon plasma at 
last scattering would free-stream to large SH multipoles 
at present, and consequently, would be easier to establish 
from CMB observations. A program of study to relate 
the BipoSH brightness fluctuations in the baryon-photon 
fluid for different physical scenarios is currently under- 
way. We have used our formalism to represent and match 
the well known case for SI violation in presence of a ho- 
mogeneous magnetic field. We illustrate how the angular 
correlations in such a case could be seeded by the dipole 
term of the generalized CMB brightness fluctuation and 
would have diagonal (£' = £) and off-diagonal (£' = £±2) 
terms. 

In summary, our work strongly motivates closer study 
of all possible SI violating phenomena and scenarios in 
the simple baryon-photon plasma, since these could po- 
tentially provide more readily observable signature of SI 
violation in the universe. It is also encouraging that 
it may have observational implications in light of the 
recent WMAP-7 discovery Q. Since, the quadrupolar 
anisotropy anomaly with non-zero BipoSH coefficients re- 
lated as A™ ~ —2A™ 21 , rule out an origin in anisotropic 
power spectrum, this may well be related to SI violations 
in the CMB brightness fluctuations. This possibility is 
also bolstered by the fact that the non-SI effect peaks 
at acoustic I ~ 200 scales pointing to some non-trivial 
physics at the last scattering surface. Extension of this 
formalism to CMB polarization should be readily possi- 
ble. The formalism and the initial conclusions are impor- 
tant and timely in light of higher precision and resolution 
CMB anisotropy and polarization data expected in near 
future, in particular, from the ongoing Planck Surveyor 
CMB mission. 



Appendix A: Bipolar coefficients for SI deviations 



The BipoSH coefficients are defined in equation (1191) . For a directional dependent primordial power spectrum as 
in equation (|16p. these Bipolar coefficients can be evaluated using the angular correlations in equation f| 1 T[) in the 
following way 

U 1 e 2 LL' 12 m MM' 

( ~\ \m' r LM r L'M' r JN r i 2 m 2 ( \T\ 

m\ra 2 m,m' 

We use the following symmetry properties of the Clebsch-Gordan coefficients in equation (8.4.10) from (48j 

r>i 2 m 2 (1 \l+li-l 2 nZ 2 -m 2 

= (- 1 )' 1 ^ 2 §;(- 1 ) rac lZ 1 [ -m ( A2 ) 
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and 

pJN / -] \t' —m' "^J plra { \o\ 

where ...t n = K^i + 1)(2^2 + l)—(2Ai + 1)]' nas been defined for convenient notational simplicity. 

We use the formula for summation of the product of four Clebsch-Gordan coefficients given in equation (8.7.26) from 

S3 

r l h £) 

EpLM nh M film piimi _ \ "> rr pl?,m 3 pLM ) tI a ill I ( \ A\ 

L -e 1 m, 1 £m U e 2 m, 2 e'm' U i'm'JN U e 2 m 2 i-m ~ 2-^1 LL UiL'i3^L' M' JN U L'M'f, 3 m 3 \ ^ *2 « ? ■ |A4j 

mim2mm' ^3^3 I ^3 t ^ J 

The Bipolar coefficients in equation (|A1[) simplifies to 

f J ^ I ] 

x E < * L M CiJ A «® A «l ■ (as) 

i i^ s {i'L'e 2 ) 1 J ^ 

{A^ ® A|^, are the Bipolar products in i and L' as defined in equation (f2"Tj) . 

For statistically isotropic primordial perturbations, the angular correlation in equation (|23[) can be written in terms 
of Bipolar coefficients as 



{A u , ) [=0 



^P (fc) £ /3 W £ A^(fc,r) (fc,r) 

«iLL' MM' 



X ( ^ Qmf'-m'C^mitaQimifm'- (A6) 

mimm' 

The summation in the above equation can be simplified using equation (8.7.17) from [HI as follows 



E/_T \m'pJN pLM pL'M' _ / -i \£±+L-\-J tt /iiM J ^ ^1 £ 

iimm' ^ 

-(-1) 1W(-1) n-t LML ,_ M , i f ^ ^ 
= (-i) < ' +i+i ' +J n ii ,(-ifC--M- { % i L [ } ■ (A7) 



L J i 



Thus the BipoSH coefficient can be simplified to 

',2 



(^) I=0 = 4. / ^Po(fc) E ft* { f / / 

^ £iLL' 

x £ A^(fc,r) [A0' r (fc,r)]*(-l) M 'C L T fi ,_ M , 

MM' 

= 4tt / ^Po(fc) E (-l) <1+£+£ ' +J n LI/i 9 w { J / ^' 

x {a^® A^,} jjv , (A8) 
where symmetry properties of the Clebsch-Gordan coefficients, equations (8.4.10) in [48| has been used. To consider 
low Bipolar deviations from statistical isotropy i.e. (L, L', J << £, £', £1), we use the asymptotic relation for Wigner-6j 
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function as given in equation (9.9.1) from [45 



L J L'\ . (-1) L +^ 



„ (-iy J+l+i 'j { e-n 



Using the above relation, the asymptotic limit to the BipoSH coefficients turn out to be 



^L(i-i^L'(£-,-£'V (^9) 



(^ J A=o « **j ^ftWEt- 1 ^^ ^ ft* ^ig^^o {A[ l£ ® a^} jjv .(aio) 



Appendix B: Generalized evolution equation for statistical isotropy breakdown 

Starting with equation (|29")1 . putting P^(fc ■ ft) = ^j^j 5Z m XfmW ^m(ft) an d evaluating the double integral using the 
equation 

<m A Y em (h) Y e±mi (h)Y e2m2 (ft)* = -]= ^-C*f ii0 Cf%£ mi , (Bl) 



the most general evolution equation becomes 

Kt^r) = ^-EH' n l ^( fcA ^) E^A^f (fc,r s )^C s , (B2) 

where 

^ _ r t 3 r t 2 r-LAf nL'W ST^ nl 3 m 3 nl 2 m 2 mo\ 

r m\m 2 m 3 m 4 m 

Using the summation formula for four Clebsch-Gordon coefficients given by equation (9.1.8) from (48[, 

E^L'M' ^ 3 D3 r iM />4mj —A A C i \<?i+W4+£'tt J ^ ^ ^ 3 I fR/H 



equation JB2| simplifies to 



ee 3 e 4 ' ^ 



E H)* n "^^£t 4(*At)A&(*,t.) C& C&° { £ * | } . (B5) 

Substituting the multipole dependent coefficients as 

C(LJJ lj2 ) = {-if H miti ^± (-l)^ +L C*f ii0 C%1, { l l } , (B6) 

the generalized evolution equation for deviations of statistical isotropy in the CMB brightness fluctuation reduces to 

&tX(k, T ) = C{L,l,h,h) Je (kAr) A^(fc,r s ) . (B7) 
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Appendix C: Diagonal terms of the asymptotic moments of CMB brightness fluctuation 



In equation (|34|) , taking e = i.e. I = R and putting £\ = R + £4 & 4 = R + 4, we can write equation (|35|) as 

, (-1) L 



4 L 4 
£ + £ 4 I t + h 



21 n. 



(CI) 



From the properties of the Clebsch- Gordon coefficient, we get 4 = £4 and hence £\ — £2- Thus in equation (|3Tj) 
substituting £ 3 = £4 — £, we get the factor 



r< ( T D D \ - ( ,V TT 2 TT TT fit x -l,l x -l < , kL r (h-£)0 r i t ( 1 ) L 

Ul{L,l,l L ) - {-l) 11^ II4 (-1) L <Mi0 -«0(^i-«)0^fn7~ 



^(.ti-tHd-t) 
" u (£i-£)(^i-^)iO 



Now using equation (8.4.10) from [48j . 



C 



eoiiO 

£0(^i-f)0 



1-^)0 _ r (ti-£)0 



ib 



(_l)«_r£j_ C (^-')o 



Using equations (8.5.34) and (8.5.42) from [Hj], we get 



4! 



~] 2 



£!(4 -4M 



{2£)\{2£4 - 2£+l)\ 
(24 + 1)! : 



and 



Thus 



2lx - 21 + 1 



(24 -21- L)!(24 -2£ + L + l) 



C4(LJ,h) 



nSn,, it 



2f 



4! 



(2^)1(2^ + 1)! (24)! (24)! 



(24 + 1)! V (24-i)! (24+i + l)! 



Substituting this in equation (|31l) . the evolution equation becomes 



(C2) 



(C3) 



(CM) 



(C5) 



(C6) 



(C7) 



where £,4 » |4 - £\. 



Appendix D: Off-diagonal terms of the asymptotic moments of CMB brightness fluctuation 
In equation (|3~4"1) , with e = i.e. £ = R 

\4 Z 4 J V2lU e ., M^-*)£(*i-40' 

From conditions of Clebsch-Gordon coefficients 

< 4 
|4 - i\ < 4 
|4-4| <£■ 

For off-diagonal terms of we consider the minimum values of 4> 4 & £ with neither equal to 0. 

4-^ = 4 
4-^ = 4 
4 - 4 = I. 
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Thus in equation (|3T1) we would get the factor 

c 2 (^,4,4) = H) £ nfn £l n £2 _, 



a ( -i\L+£ 1 +t 2 



c 



{lx-l){l^t) 



Now using the symmetries of Clebsch- Gordon coefficients given by equation (8.5.34) from [48|, we get 



(D2) 



r {lx-C)Q _ r {ti-l)0 



Similarly 



C 



e 2 o 

eo(t 2 -e)a 



(-l)^i! 
t\{ti-l)\ 



(2l)!(24-2l+l) 
(24 + 1)! 



11^2 



n <2 4! 
^!(4 - 1)\ 



(2£)!(2£ 2 -2^ + l)! 



(24 + 1) 



Using equation (8.5.37) from |48[ we get 



^{l 2 -l)U 2 -l)LL ~ °> 



8i 2 +L,e-L- 



Thus 



l + *2 



(D3) 



(D4) 



(2^!^!£ 2 ! 



(24 - 2^ + l)!(2^ - 2^+1)! 



(2^ + 1)1(2^ + 1)! 



n 



£«i£ 2 



V2I (ei)^e 1 -ey.(£ 2 -£)\ 

Substituting this in equation (|31l) . the evolution equation becomes 

Ai#,(fc,r) = ^C 2 (^,4,4) ii(fcAr) A££ Va _,(fc,r.) fc,*. 



4+4-2^ + 1 



n^_ f V 4 + 4 + 1 



(D5) 



(D6) 



(D7) 



where I, £ n » \t n - £\, L with n = 1, 2. 
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